An introduction to the spectrum, symmetries, and dynamics of Heisenberg spins- 1/2 chains 
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Quantum spin chains are prototype quantum many-body systems. They are employed in the description of 
various complex physical phenomena. The goal of this paper is to provide an introduction to the subject by 
focusing on the time evolution of a Heisenberg spin- 1/2 chain and interpreting the results based on the analysis 
of the eigenvalues, eigenstates, and symmetries of the system. We make available online all computer codes 
used to obtain our data. 



I. INTRODUCTION 

The fascination with magnets can be traced back as far as 
the antiquity in China, but it was only with the discovery of 
spins [[0, 12D that we developed a better understanding of mag- 
netism. Insulating solids with magnetic properties, in partic- 
ular, can be viewed as lattices of atomic or ionic magnetic 
moments, each localized to one site JH-Hl]. The total magnetic 
moment of the atom or ion depends, in general, on the spins 
of the electrons in incomplete shells and their orbital angular 
momenta. Here we refer to this total angular momenta simply 
as spin. The interactions between the spins may lead to col- 
lective behaviors with macroscopic effects, such as ferromag- 
netism, where the spins line up parallel to each other, and an- 
tiferromagnetism, where neighboring spins point in opposite 
directions. The source of such spontaneous magnetization is 
the so-called exchange interaction, introduced by Heisenberg 
and Dirac in the end of the 1920's JH0]. This interaction has 
a quantum mechanical origin. It is the manifestation of the 
Coulomb repulsion between the electrons and the Pauli ex- 
clusion principle, being therefore strong and short range. A 
magnetic dipolar interaction is also present, but it is too small 
to explain magnetism at room temperature. 

A general form of the exchange interaction is given by the 
Heisenberg model. This is one of the most important models 
in magnetism iH-HI and also in studies of exactly solved quan- 
tum many-body systems J^] . Despite being a simplified theo- 
retical model, it describes quantitatively well some magnetic 
compounds, such as copper oxides Q. In some of these real 
systems, especially the quasi-one-dimensional ones, unusual 
high transport of heat has been verified 0$^, which further 
motivated the interest in the model. Such anomalous trans- 
port behavior has been associated with a macroscopic number 
of conserved quantities characterizing the one-dimensional 
Heisenberg spins- 1/2 model with coupling between nearest- 
neighbor sites only lioll . In fact,under these conditions, the 
model is integrable and can be analytically solved with the 
Bethe ansatz method lfTTl[T2ll . 

The Heisenberg model finds applications in several other 
contexts. It is a key model in studies of quantum phase transi- 
tion ||l3ll . of superconductivity ft^, as well as of the dynamics 
and thermalization in correlated one-dimensional lattice sys- 
tems Il4tl . In quantum information, Heisenberg systems are 
used as models for quantum computers, each spin- 1/2 rep- 
resenting a quantum bit (qubit), in the analysis of entangle- 
ment IU5I1 . and in methods to transfer information in a con- 
trollable way HEl. In the presence of impurities, disorder, 



or couplings beyond nearest-neighbors, the system becomes 
non-integrable and has been employed in the characterization 
of the crossover from integrability to quantum chaos lfl8l - l2Tll . 
Furthermore, there has been atternpts to simulate this model 
with cold gases in optical lattices 112211 . 

Optical lattices are crystals of light. Laser beams propagat- 
ing in opposite directions result in standing waves that confine 
ultracold atoms to small regions, the atoms playing the role 
of electrons in solid crystals I23I1 . These systems are highly 
controllable, which allows for the simulation of condensed 
matter models not easily accessible with real solid state sys- 
tems. Moreover, they are weakly coupled to the environment, 
which makes it possible to study their evolution for a long 
time. These factors combined make optical lattices essential 
tools to advance of our understanding of quantum many-body 
systems far from equilibrium. The behavior of nonequilib- 
rium systems is an outstanding challenge at the forefront of 
physics. 

Motivated by the widely spread interest in spin systems and 
in the out-of-equilibrium properties of quantum many-body, 
we study here the factors that may limit the time evolution 
of the one-dimensional Heisenberg spin- 1/2 model with cou- 
plings between nearest-neighbor sites only. We focus on the 
effects of an anisotropy parameter and on the symmetries of 
the system. We rely on the analysis of the eigenvalues and 
eigenstates of the system to anticipate its dynamics. Our pre- 
dictions are then confirmed with actual numerical results for 
the time evolution. Since our studies require all eigenvalues 
and eigenstates, we use exact diagonalization. 

The paper is organized as follows. Section II provides a 
detailed description of the Hamiltonian of the system. Section 
III analyzes the diagonal elements of the Hamiltonian, as well 
as its eigenvalues and eigenstates. Section IV investigates the 
time evolution of the system. Its symmetries and how they can 
constrain the system dynamics are discussed in Sec. V. 



II. SPIN-1/2 CHAIN 

The description of a system composed of a single spin- 
1/2 requires the use of the spin operators S x,v,z = a x ' v ' z /2, 
where 
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are the Pauli matrices and h has been set to 1 . The quantum 
state of the spin is represented by a two-component vector, 



2 



known as the spinor. This state is commonly written in terms 
of basis vectors corresponding to the eigenstates of S z , which 
are the spin pointing up in the z direction | f) = ( Q ) and the 
down-spin | \) = (.) . The eigenvalue associated with | f) 
is +1/2 and that of | J,) is -1/2, which justifies referring to the 
first as the excitation. 

Here, we study a one-dimensional system (chain) com- 
posed of L coupled spins- 1/2 described by the Heisenberg 
model, 



H = 



L-l 

E 

Tl=l 



J 1 S~, S t 



qy qv 



(1) 



The operators S^' v ' z act only on the spin placed on site n. 
The couplings are limited to nearest-neighbor spins; J is 
the strength of the flip-flop term S%S% +1 + J z is 

the strength of the Ising interaction S Z S Z +1 , and the ratio 
A = J z j J is the anisotropy parameter. The model is isotropic 
when A = 1, in which case it is known as the XXX model, 
and it is anisotropic when A ^ 1, usually referred to as the 
XXZ model (XYZ also exists when the coupling strengths in 
the three directions are different). A natural basis for the sys- 
tem is the set of 2 L states where the spin on each site is either 
pointing up or down, such as | iiT2T3 • ■ • II)- In quantum in- 
formation, these states are known as quantum computational 
basis vectors. We refer to them as site-basis vectors. 

The flip-flop term interchanges the position of neighboring 
up and down spins according to 



J[S^S, 



n+l 



^n^n+l) I tnln 



+ 1/ 



(J/2)| i 



n I n+l I 



It therefore couples site-basis vectors that differ only by the 
orientation of the spins in two adjacent sites. In this basis, it 
constitutes the off-diagonal elements of the Hamiltonian ma- 
trix. This term has a key role in the evolution of the system by 
moving the excitations through the chain. 

In the case of open boundary conditions (open chain), as in 
Eq. ([TJ, where the sum goes from site n = 1 to site L — 1, 
an excitation on site 1 can move only to site 2 and from site 
L to L — 1. The scenario of a ring (closed chain), where an 
excitation on site L can also move to site 1 corresponds to 
closed boundary conditions and will be discussed briefly here. 

The Ising interaction contributes to the diagonal part of the 
Hamiltonian matrix written in the site-basis. It causes a pair 
of adjacent parallel spins to have different energy from a pair 
of anti-parallel spins, because 



JzSnSn+l\ tntn+l) — + (Jz/4)\ tntn+l/j 



while 



JzS„Sj l+1 \ tnin+l) — — (J*/4)| tnlri+l)- 



(2) 
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A Hamiltonian containing only the Ising interaction, H zz = 
J Z S^S^ +1 , constitutes the Ising model and was employed in 
the first attempts to describe the phase transition from para- 
magnetism to ferromagnetism [24]. As one can infer from 
Eqs. (O and (01, the ground state of this model depends on 



the sign of the interaction strength; it is ferromagnetic, with 
all spins aligned in the same direction, when J z < 0, and 
it shows an antiferromagnetic arrangement with antiparallel 
neighboring spins when J z > 0. 

The state in which all spins align in the same direction is 
also an eigenstate of the Heisenberg model, because the flip- 
flop term has no effect on it. When the Heisenberg model 
is ferromagnetic ( J z < 0), this state is a ground state. For 
the antiferromagnetic Heisenberg model (J z > 0), on the 
other hand, the ground state is more complicated than the sim- 
ple configuration of antiparallel spins (details about it may be 
found, for instance, in Refs. J5I). 

Here, we are interested not only in the ground state, but 
in all eigenvalues and eigenstates of the finite antiferromag- 
netic XXZ model with open boundary conditions. They are 
computed numerically and used to compare the static and dy- 
namic properties of the system in two scenarios, when A < 1 
and when A>1. We address the role of the anisotropy pa- 
rameter, as well as border effects and symmetries. 

Our Hamiltonian commutes with the total spin in the z di- 
rection, S z = ^2n=i^n' tr,at i s ' [H,S Z ] = 0. This means 
that the system is invariant by a rotation around the z-axis, or 
equivalently, it conserves S z . As a result, the Hamiltonian ma- 
trix of a system with L sites is composed of L + 1 independent 
blocks (or subspaces), each with a fixed number N € [0, L] 
of up-spins. Therefore, even though the total dimension of the 
Hilbert space is 2 L , we can diagonalize a single subspace at 
a time, each of dimension D = (^) . Other symmetries may 
also be present, but this discussion is left for Sec. V. In the 
next two sections, we focus on the analysis of the spectrum of 
a particular S z -subspace and the level of derealization of its 
eigenstates, as a way to predict the dynamics of the system, 
and then confirm our expectations by investigating the actual 
time evolution of different initial states. 



III. SPECTRUM 

Before diagonalizing the XXZ Hamiltonian matrix, which 
we write in the site-basis, let us first look at its diagonal ele- 
ments. They form bands of degenerate energies. These bands 
are determined by the number of pairs of adjacent parallel 
spins in the basis vectors. For each S z -subspace, the larger 
the number of pairs, the larger the energy of the basis vector. 
For example, in an open chain with L = 4 and N = 2, the 
highest energy, J A/4, occurs for the states with two pairs of 
parallel spins, | ttii) and I lift) ■ The band that precedes 
this one in energy has the states with only one pair of paral- 
lel spins, I tilt) and I lttl)> yielding an energy of - J A/4, 
while the states of the band with the lowest energy, — 3 JA/4, 
have no pairs of parallel spins, | tltt) an d I tttt)- One 
sees that the energy difference between consecutive bands is 
J A /2, since we move down in energy by breaking a pair, thus 
adding the factor JA/4 one less time and subtracting it one 
more time. In an open chain, where there are L — 1 coupling 
bonds, the general expression for the energy of each band is 
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therefore 



Br = [2p-(L-l)]JA/4, 



(4) 



where p is the total number of pairs of adjacent parallel spins. 

In a closed chain, on the other hand, the energy difference 
between successive bands is J A. In this case, there are L 
bonds and always an even number of antiparallel pairs, be- 
cause there is no border to absorb any of them. We move 
down in energy by breaking necessarily two pairs of parallel 
spins, the factor J A/4 being thus added two less times and 
subtracted two more times. The diagonal energies are then 
given by 



E 



closed 



[2p-L}JA/4. 



(5) 



Clearly, the closed chain has fewer bands than the open one, 
as shown in Table Q] (Note: [x J is the integer part of x). 



TABLE I: Number of energy bands formed with the diagonal ele- 
ments of the XXZ Hamiltonian written in the site-basis. 
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FIG. 1 : (Color online) Histograms of the diagonal elements [panels 
(a) and (b)] and of the eigenvalues [panels (c) and (d)], and IPR aver- 
aged over all eigenstates vs A [panel (e)], for the antiferromagnetic 
XXZ Hamiltonian with open boundaries, L — 10, and N = 5. The 
site-basis is used. Panels (a) and (c): A = 0.5; panels (b) and (d): 
A = 10. Panel (c): bin width = 0.2 and panel (d): bin width = 1.0. 
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The bands for the diagonal energies of an open chain with 
an even number of sites and N = L/2 are illustrated with 
histograms in Figs. Q] (a) and (b) for two values of A. The 
energies range from — {L — 1) JA/4 to (L — 3)JA/4. The 
L — 1 bands are symmetric, so the number of states tjb in 
band B is equal to that in band L — B, that is j]b = i]l-b- 
The number of states in each band grows as we approach the 
middle of the spectrum and is given by 



VB = tjl-b 



211 



fc=i 



N - [k/2\ 



NS, 



fe.i 



(6) 



where B < L/2. The least populated bands are the ones in 
the extremes containing only two states each, r/i = t)l-\ = 
2. If L mod 4 7^ 0, as in the figure, B = L/2 is the most 
populated band, whereas if L is divisible by 4, the three bands 
in the middle are the most populated with rj L / 2 -i = Vl/2 = 
Vl /2+1 ■ (Illustrations for the bands of the closed chain can be 
obtained with the codes in |[25lo . 

After diagonalization, the structure of the spectrum depends 
on the interplay between the Ising interaction and the flip-flop 
term. FiguresQ~](c) and (d) show the histograms for the eigen- 
values of the Hamiltonian (fl} with the same parameters con- 
sidered in the top panels. In Fig.[T](c), where A < 1, the band 
structure is lost. This happens because the energy difference 
between the basis vectors is < J, so the flip-flop term can cou- 
ple intra- and also inter-band states, broadening significantly 
the range of energy values. In contrast, the band structure is 
preserved in Fig.[TJ(d). There, since A ^> 1, states from dif- 
ferent bands are too far off-resonance and the flip-flop term 



can effectively couple only states belonging to the same band. 
Each band then acquires a small width, which does not erase 
the energy gap between them ll26ll . 

The competition between the Ising and the flip-flop term of 
the Hamiltonian is, of course, reflected also in the structure of 
the eigenstates. As the Ising interaction increases, limiting the 
role of the flip-flop term, the eigenstates become less spread in 
the site-basis. This can be quantified, for example, with the so- 
called inverse participation ratio (IPR) IT271 l28ll . Consider an 
eigenstate written in the basis vectors \(f>k) as = 



<.'> 



\4>k)- IPR is defined as 

jU) = 



1 



yD 



(7) 



This quantity is proportional to the number of basis vectors 
which contribute to each eigenstate. It is small when the state 
is localized and large when the state is delocalized in the cho- 
sen basis. 

Figure Q](e) shows IPR averaged over all eigenstates, (I), 
for various values of the anisotropy parameter. The maximum 
derealization occurs at A = 0. As the anisotropy increases, 
(I) decays monotonically until the energy bands cease over- 
lapping and (7) approaches a constant value ll29ll . 



IV. DYNAMICS 

We now analyze the time evolution of different initial states, 
each corresponding to a specific site-basis vector, ^(O)) = 
\4>k)- The source of the dynamics is the flip-flop term, which 
couples \4>k) with other states, transforming into an 

evolving superposition of site-basis vectors. From the results 
for the eigenvalues and eigenstates described in the previous 
section, we expect the initial state to spread over several basis 
vectors when A < 1, whereas the dynamics should be con- 
fined to states belonging to the same energy band as ^(O)) 
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when A ^> 1. To confirm this prediction, we study two quan- 
tities, the magnetization of each site, 



M n (t) = (*(t)|£*|*(f)>, 



(8) 



and the probability Pi(t) for finding a basis vector \<j){) at in- 
stants of times t. Since 

i=i 

= EfE4 j) «pv^ii^ 
i=i \i=i 

the probability is given by 



fi(*) 



D 



i=i 



(9) 



Note that the coefficients a's are real, because the Hamiltonian 
(dJ is a real symmetric matrix. 

The top panels of Fig.[2]show the magnetization of each site 
of an open chain that has a single excitation initially placed on 
site 1, |*(0)) = | t-W- • • ■ D- In P anel ( a ). where A = 0.5, 
we see the up-spin leaving the edge and gradually spreading 
through the chain by hopping successively from one site to the 
next in intervals of time ~ J" 1 . The probability of finding it 
on a single site decreases from site 2 to L — 1. The excitation 
then finally reaches the other edge with high probability. The 
preference for the edges is a border effect that appears only 
in an open chain and is caused by the resonance between the 
two border states | t44- ■ • ■ i) and I ii ■ ■ ■ IT)- Despite this 
resonance, the other states also take part in the dynamics, be- 
cause their energy differ from that of the border states by only 
JA/2 = 0.25J. For comparison, we show in panel (b) the 
case where this energy difference is large, JA/2 = 5 J. The 
intermediate states between the edges are now virtual and the 
spins on sites 2 < n < L—l remain pointing down throughout 
the evolution. Only the border states are effectively coupled 
and in a high order of perturbation theory, which explains the 
long time it takes for the excitation to move from n = 1 to 
n = L. Clearly, this time becomes even longer as the size of 
the chain increases. 

The different time scales associated with the order of per- 
turbation theory in which the states are effectively coupled are 
well illustrated in Fig.|2](c). There we have 6 sites and 2 exci- 
tations initially placed away from the borders. The initial state 
1^(0)) = I ITtm) is m resonance with the basis vectors 
I -UTt-U). I i-Utti)- and | t-UUt)- These four states be- 
long to the same energy band, which remains separated from 
the other bands once the flip-flop term is turned on, because 
the anisotropy considered is large, A = 10. Thus, we expect 
only these four states to be able to mix. 

State | mm) couples with state | mm) in second or- 
der of perturbation theory via an intermediate transition where 
the pair splits into the virtual state | ^tltll) ar, d then recom- 




Jt(xlO ) 



FIG. 2: (Color online) Top panels: Magnetization of each site ver- 
sus time; initial state | tlllll); A = 0.5 [panel (a)] and A = 10 
[panel (b)]. Bottom panels: probability in time to find a specific basis 
vector, A = 10. Panel (c): Initial state | 4TT444-) (solid line) cou- 
ples with | 44tt44) (circles) , | 444-tM-) (triangles), and | t-U-Ut) 
(plus). Panel (d): initial state | ttt-Ui) couples only with | 444ttt)- 
All panels: open chain, L — 6. 



bines again. The effective coupling strength obtained by per- 
turbation theory is therefore 



eff ' 



mtm \h, v \ mm> (mm \h xv \ mm) _ j 



I4.n-W.4-> 



E 



41 il-i-L) 



where H xy = J(S^Sf l+1 + S^S^ +1 ). The same coupling 
strength is found between states | mm) and | mm)- 
These states then hybridize at t ~ J^- 1 = A/ J. In contrast, 
the effective coupling between | mi-U) and | mm) oc- 
curs in fourth order of perturbation theory. As a result, it takes 
much longer for the latter state to take part in the dynamics, 
as clearly seen in Fig.[2](c). 

The combination of large anisotropy and border effects can 
slow down the evolution of the system significantly to the 
point that an initial state may look stationary IB0l - l32ll . This 
is the case shown in Fig. |2] (d), where A = 10 and the initial 
state is |*(0)) = | tttm)- This state can hybridize only 
with | 444-TtT)- The communication between the two states 
occurs in a very high order of perturbation theory, the exci- 
tations taking a long time to cross the chain from the left to 
the right side. This makes the initial state look frozen for a 
long time. Notice that in a closed chain, where the borders are 
absent, the dynamics would be faster. In this case, the bound 
triple of excitations Ttf would not be restricted to the edges, 
but would move together through the whole chain. 

Separated energy bands caused by large anisotropy and the 
presence of borders can then limit the dynamics of the system 
to a portion of the Hilbert space. Another restrictive factor is 
the symmetries of the system, as we discuss next. 



V. SYMMETRIES 

Operators that commute with the Hamiltonian have two im- 
portant properties: their eigenstates are also eigenstates of 
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H B33I1 and they represent physical quantities that are con- 
served i34ll . The latter property comes from the fact that 
the expectation value of a conserved quantity, (O), does not 
change in time, so 



d(Q) 
dt 



-{HO - OH) = =>- [H, O] = 0, 



The way to find such constants of motion is by looking for 
the symmetries of the system. According to Noether's theo- 
rem 0511 . the invariance of the Hamiltonian under a symmetry 
operation must necessarily lead to a conserved quantity. For 
example, invariance of H under translation in space leads to 
conservation of linear momentum, and invariance of H under 
translation in time leads to conservation of total energy. 

We have already encountered a conserved quantity of our 
system, the total spin in the z-direction. As we saw when 
studying the dynamics, the value of S z of the initial state is 
conserved throughout the evolution. If the initial state has N 
up-spins, then all the states that take part in its evolution must 
have the same number of up-spins. But our system shows 
additional symmetries, as we describe next. 

Hamiltonian dTJ is invariant under reflection, which leads to 
conservation of parity, that is, H commutes with an operator 
II that transforms the spin operators as S n — > Sl-u- This 
may be better understood by imagining a mirror at one edge 
of the chain. If parity is conserved, the probability of each 
basis vector in the eigenstate is equal to that of its reflection. 
For example, suppose we have L = 4 and one excitation. The 
eigenstates are given by 



tm>- 



Itll)- 



mi)- 



,0')i 



lilt) 



and the probability amplitudes are either a\ 



U) 



and 



4° 



,U) 

3 

,00 



for even parity, IT 



-1, or a\ 



00 



-a( and 



— a 3 for odd parity, IT = — 1. 



If L is even and N = L/2, our Hamiltonian is also invariant 
under a global tt rotation around the x axis. We denote the 
operator that realizes this rotation by R%. As an example, 
suppose we have L = 4 and N = 2. The eigenstate 
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TABLE II: The first column numerates the eigenstates \ipj) of an 
open XXZ chain with A = 0.4, L = 6, and N = 3. The sec- 
ond and third columns give, respectively, the eigenvalues of IT and 
R% of these eigenstates. Fourth, fifth, sixth, and seventh columns: 
probability amplitudes a,j of the initial states {ToJ written as super- 
positions of the eigenstates of the Hamiltonian, |^E'a,s,c,d(0)) = 
a/j ' B ' C ' D '\tpj). The eigenvalues offlandit^ for these initial states 
are shown in the second and third rows. 
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-1. 



, in which 
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There are two other symmetries, which we will not explore 
here. One is conservation of total spin, St = J2 n ^ n > wn i c h 
occurs only in isotropic systems, where [H,S T ] = 0. The 
other is conservation of momentum, which happens in the 
closed chain due to its invariance by a translation in space. 

In this work, we focus on IT and R* and analyze how they 
affect the dynamics of the system. For this, we consider four 



different initial states, 



I*a(0)> = 


(1 ItttU) 


+ llltttl))/^ 




I*b(0)) = 


(1 ItttU) 


-llltttl))/^ 


Fi 


1*0(0)} = 


(1 ItttU) 


+ ltUltt))/^ 


Fi 


\*d(0)) = 


(1 ItttU) 


- 1 mm) + 






1 tUItt) 


-ttUlt))/2. 





(10) 



The first two states are not eigenstates of R%, but parity is 
well defined, ^^(O)) has even parity and |^s(0)) has odd 
parity. State j'T'c^O)), on the other hand, is an eigenstate of the 
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operator R* with eigenvalue +1, but not of II. The last state, 
|5 , d(0)) is an eigenstate of both II and R%, with eigenvalues 
— 1 and +1, respectively. 

In Table |nl we write each initial state as a linear superposi- 
tion of the eigenstates of the system, so that we can investigate 
which eigenstates can take part in the evolution. The first col- 
umn numerates the eigenstates. The second and third columns 
give the eigenvalues of IT and R%, respectively, for each eigen- 
state. The fourth, fifth, sixth, and seventh columns show the 
values of the probability amplitudes of the eigenstates for the 
initial states |*a(0)), |*b(0)>, |*c(0)), and |#d(0)>, re- 
spectively. The eigenvalues of II and for each initial state 
are shown in the second and third rows. 

The table makes it clear that only the eigenstates with the 
same symmetries as the initial state can contribute to the evo- 
lution of the latter. For|\&yi(0)) and ^^(O)), eigenstates with 
both values of R* are seen, but parity is strictly conserved. In 
the first case, the probability amplitudes of all odd eigenstates 
are zero and for the second state, a,j = for all eigenstates 
with II = +1. For state \fyc(ty}> eigenstates with both pari- 
ties are part of the superposition, but aj ^ only for eigen- 
states with R% = +1. The last state, |5 , d(0)), has both sym- 
metries, so only three of the 20 eigenstates have ctj ^ 0, those 
with II = — 1 and R* = +1 simultaneously. These three 
eigenstates, ^3,^5, and ipie, are more localized in the site- 
basis than all other eigenstates. This is because superpositions 
of some of the basis-vectors cannot satisfy both II = — 1 and 
R* = +1 at the same time. For example, the effect of the 

two operators IT and R% on bi\ ttt-UI) + M 444-1TT) is m e 
same, so they cannot give different eigenvalues. These two 
basis vectors (and similarly for others) cannot, therefore, be 



part of the contributing eigenstates -03, and V>i6- 

In summary, the larger the number of conserved quantities, 
the smaller the invariant subspaces. The time evolution of an 
initial state with many constants of motion is therefore more 
constrained in the Hilbert space. 



VI. DISCUSSION 

We make available online, in l25ll . all computer codes used 
to obtain our data, as well as detailed explanations about them. 
Students and professors should have no difficulty to reproduce 
our results and to explore further questions. The studies de- 
scribed here can constitute an entire summer project, as in our 
case, can set the basis for a senior thesis, or give ideas for as- 
signments in courses about Quantum Mechanics. By the time 
of completion of this work, the two first authors of this pa- 
per were undergraduate students who had not had a course on 
Quantum Mechanics, but had a solid background on Linear 
Algebra. 
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